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Abstract 

In usual scale-free networks of Barabasi-Albert type, a newly added node selects 
randomly m neighbors from the already existing network nodes, proportionally 
to the number of links these had before. Then the number N(k) of nodes with k 
links each decays as 1/fc 7 where 7 = 3 is universal, i.e. independent of m. Now 
we use a limited directedness in the construction of the network, as a result of 
which the exponent 7 decreases from 3 to 2 for increasing m. 
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Introduction: 

The Barabasi-Albert network is growing such that the probability of a new site 
to be connected to one of the already existing sites, is proportional to the number 
of previous connections to this already existing site: The rich gets richer, in this 
way, each new site selects exactly m old sites as neighbors. 

In directed Barabasi-Albert networks [1], the network itself was produced in 
the standard way, but then the neighbor relations were such that if A has B as 
a neighbor, B in general does not have A as a neighbor [2]. 

The undirected Barabasi-Albert network [3], usually is grown in the same 
way, but then the neighbor relations were such that if A has B as a neighbor, 
B has A as a neighbor [2]. 

Now we try to introduce a directedness already in the construction of the 
network[3]. In the undirected Barabasi-Albert network [4-5], if a new node 
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selects m old nodes as neighbor, then the m old nodes are added to the Kertesz 
list, and the new node is also added m times to that list. 

Connections are made with m randomly selected elements of that list. If 
one would only add the old nodes to the list, then only the initial core can be 
selected as neighbors, which is not interesting. But if one adds the m old nodes 
plus only once (and not m times) the new node, then one has a semi-directed 
network. In this modified Barabasi- Albert model one can (but we don't do that 
here) put in neighbor relations which are directed or which are undirected. 

Data and Simulation: 

We use the FORTRAN program[6] , as in appendix, with different m= 1, 2, 3, . . . , 10, 
and the maxtime = 41,000,000 nodes. 

Distribution of N(K) of nodes with K links 
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Figure 1: n(k) versus k links at number of neighbors m = 1,2, and 10, from 
left to right, maxtime=40000000 in log- log plot. 

Figure(l) shows the number n(k) of nodes with k links each. Changes are 
seen with increasing number of neighbors m, and the curvature shows that the 
exponent should be determined from the region on large k only. 

We are going to plot the slopes versus 1/m to get the second plot , which 
makes clearer the possible extrapolation towards infinite m(m = oo, 1/m = 0) 
as shown in figure 2. 
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Figure 2: Slopes versus 1/m for number of neighbors m = 1, 2, 3, . . . , 10, max- 
time=41000000. 

Conclusions: 

The exponent gamma for the decay of the number n(k) of nodes with k links 
each, changes with increasing number of neighbors m, the exponent should be 
determined from the region on large k only. 

Fig. 2 shows 7 — 1 since we summed all k into bins limited by powers of two, 
e.g, 4, 5, 6, 7 is one of the many points onto which we fitted the exponent 7 — 1. 
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Appendix: 



parameter (nrun=l ,maxtime=40000000 ,m=2 , iseed=l ,max=maxtime+m, 
1 length=l+(l+m) *maxtime+m* (m-1) ) 

integer*8 ibm 
real*8 factor 

dimension k (max) , nk(131070) , list (length) , nklog(30) 
data nk/131070*0/,nklog/30*0/ 
OPEN (UNIT=60,FILE= , N(K)M2.DAT') 

WRITE (60,*) '#(nrun, maxtime, m, iseed) ' ,nrun, maxtime, m,iseed 
ibm=2*iseed-l 

f actor=(0 . 25/2147483648 . OdO) /2147483648 . OdO 
fac=l. 0/0. 69314 
c factor=0. 5/2147483648. OdO 

do 5 irun=l,nmn 
do 3 i=l,m 
do 7 j=(i-l)*(m-l)+l, (i-l)*(m-l)+m-l 
7 list(j)=i 

3 k(i)=m-l 
L=m* (m-1) 
if(m.eq.l) then 

L=l 

list(l)=l 
endif 

c All m initial sites are connected with each other 

do 1 n=m+l,max 
do 2 new=l,m 

4 ibm=ibm*16807 
j=l+(ibm*factor+0.5)*L 

if ( j . le . . or . j . gt .L) goto 4 
j=list(j) 
list (L+new)=j 
2 k(j)=k(j)+l 
list(L+m+l)=n 
L=L+m+l 
1 k(n)=l 

WRITE (60,*), '# (irun )', irun 
C print *,'# (irun )', irun 

do 5 i=l,max 
k(i)=min0(k(i) ,131070) 

5 nk(k(i))=nk(k(i))+l 
do 6 i=l, 131070 

6 if (nk(i) .gt.0) write (60,*), i,nk(i) 
do 9 i=l, 131070 

j=alog(f loat(i))*f ac 

9 nklog(j)=nklog(j)+nk(i) 
do 10 j=0,18 

10 WRITE (60,*), sqrt(2.0)*2**j ,nklog(j) , j 
stop 

end 
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